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Abstract: We study the time dependence of the entanglement entropy of disjoint in¬ 
tervals following a global quantum quench in (l-l-l)-dimensional CFTs at large-c with a 
sparse spectrum. The result agrees with a holographic calculation but differs from the 
free held theory answer. In particular, a simple model of free quasiparticle propagation 
is not adequate for CFTs with a holographic dual. We elaborate on the entangle¬ 
ment tsunami proposal of Liu and Suh and show how it can be used to reproduce the 
holographic answer. 
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1 Introduction 

Computation of entanglement entropy in quantum held theory is a topic of growing 
interest. Unfortunately, in many situations it is an extremely difficult quantity to 
compute. For strongly coupled systems, one of our best analytical tools is holography. 
In a holographic theory, entanglement entropy is evaluated as the area of an extremal 
surface in the dual geometry [1, 2]. In some cases, especially in (l+l)-dimensionaI 
CFTs, purely held-theoretical arguments can be given for the form of the entanglement 
entropy. 

We will be using entanglement entropy to study the approach to equilibrium fol¬ 
lowing a global quantum quench. Calabrese and Cardy [3, 4] famously showed that, in 
(l-l-l)-dimensions, the entanglement entropy of an interval increases linearly with time 
until it saturates at the thermal value. The held-theoretical results could be reproduced 
by a simple intuitive model: the growth and saturation of the entanglement entropy 
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was effectively modeled by the free propagation of quasiparticle excitations, which we 
will review below. 

This linear growth in (l+l)-dimensions was conhrmed in holographic calculations, 
hrst numerically [5], and then analytically [6, 7], where it was also extended to higher 
dimensions (see also [8]). In [9, 10], it was noted that a free-streaming quasiparticle 
picture was inadequate to explain the rate of growth beyond the (l-l-l)-dimensional 
case and it was suggested that the linear growth could be thought of in terms of 
an “entanglement tsunami,” represented by an effective wavefront which propagates 
into the region under consideration and entangles interior degrees of freedom with 
exterior degrees of freedom. An underlying free-streaming quasiparticle model for the 
entanglement tsunami yielded a wavefront velocity smaller than the velocity calculated 
holographically in dimensions higher than (H-l). In (1-1-1) dimensions, the holographic 
calculation and quasiparticle model both said that the effective wavefront moved at 
the speed of light. At this level, then, it seemed that quasiparticles were sufficient to 
explain the time dependence of the entanglement entropy of a single interval in (H-l) 
dimensions, even for holographic CFTs which are strongly interacting. 

However, when multiple disjoint intervals are considered, there is a qualitative 
difference between the holographic calculation and the quasiparticle model of Calabrese 
and Cardy (as noted in [11, 12], and also present in [13]): the holographic calculation 
gives a non-decreasing entropy, while the quasiparticle model has both increasing and 
decreasing behavior as a function of time. Therefore already in (l-l-l)-dimensions there 
is a need to replace the quasiparticle picture. We will elaborate on the entanglement 
tsunami proposal, showing how it can be used as a rule to calculate the entanglement 
entropy of one or two intervals in holographic CFTs and provide a natural upper bound 
for more than two intervals. We will make no attempt to derive the entanglement 
tsunami from underlying physical excitations, quasiparticle-like or otherwise. 

The remainder of this paper is organized as follows. In Section 2 we will review 
the setup for the problem and the results from both the quasiparticle and holographic 
viewpoints. In Section 3 we will provide a CFT argument for the holographic result, 
showing how CFTs with large-c and a sparse spectrum differ from weakly coupled CFTs, 
and in particular do not follow the quasiparticle prediction. In Section 4 we will propose 
an entanglement tsunami prescription for the entanglement entropy as a function of 
time. Finally, in Section 5 we will describe how these results may be extended to higher 
dimensions, where there are still some unresolved issues, and speculate on possible 
derivations of the entanglement tsunami from interactions. 
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2 Setup and Review 


2.1 Quasiparticle Model 

We are interested in the entanglement entropy of a subregion A of the real line following 
a global quench in a CFT. Note that A is always the union of a collection of disjoint 
intervals. Let A'^ denote the complement of A. Globally, the system is in a pure state 
Id'), and the entanglement entropy for the region A as given by 

S'=-Tr Pa log Pa, with pa = Wa^^ |d')(d'|. (2.1) 


|d') is a time-dependent state, and so S will be time-dependent. A global quench is 
dehned by beginning in the vacuum state of one theory, and then suddenly changing 
the Hamiltonian to that of a different theory. The result is that the initial state from 
the point of view of the second theory (the CFT) is a highly excited state, but has a 
simple entanglement structure. We are interested in the time dependence of the hnite 
part of the entanglement entropy, so we will assume that the UV-divergent parts can 
be subtracted in a consistent way, whether or not they are modihed by the quench. 
Henceforth when we refer to the entropy we will mean the hnite part of the entropy, or 
the vacuum-subtracted entropy. At t = 0, then, we have S = 0. 

At late times, the system will effectively thermalize and we should hnd S'(t —?• 
oo) = ‘S'therm = -SeqVol(A). In (l-|-l)-dimensions Vol(A) is just the sum of the lengths 
of the intervals that make up A, but we are emphasizing that the thermal entropy is 
extensive with the volume of the system. Seq is the thermal entropy density, which is a 
property of the state. At intermediate times, either a CFT calculation or a holographic 
calculation can be used to describe the transition from zero entropy to the thermal 
result. We will discuss both of these calculations below, but for now we will record the 
expected answer as predicted by the quasiparticle model. 

Calabrese and Cardy showed [3, 4] that, following a global quench, the time- 
dependence of the entropy can be effectively modeled by the propagation of entangled 
quasiparticles, at least for weakly coupled CFTs. At the time of the quench, we imagine 
that a uniform density of EPR pairs of quasiparticles are produced, where each pair 
begins localized at a point and consists of a left-mover and right-mover. These quasi¬ 
particles move in opposite directions at the speed of light, and there are no interactions 
between pairs. To compute the entanglement entropy of A, we only have to count the 
number of unpaired particles in the region at any given time: 


S{t) (X / dx' / dx'‘ 
Jx'eA Jx"eA‘^ 


dx < 6{x' — x — t)6{x''—x+t) + S{x' — x+t)6{x"-x- 


( 2 . 2 ) 


3 


Quasiparticle prediction Holographic prediction 



Figure 1. The entropy production as a function of time for a region consisting of two disjoint 
intervals of length L, separated by a distance R > L. The quasiparticle model (left) shows 
decreasing behavior between 2t = R and 2t = L -\- R. The holographic calculation (right) 
is monotonically increasing before saturation at 2t = L, after which the entropy remains 
constant. 


The constant of proportionality is related to the initial density of EPR pairs, which 
determines Seq. When A consists of a single interval of length L, computing the above 
integral gives the following time dependence: 


S{t) = 2Seq X 




t < - 


L 
2 ’ 



(2.3) 


First there is a linear growth phase, and then saturation at the thermal value. 

The behavior is a little more complicated when A consists of two disjoint intervals. 
For simplicity, consider the case where the two intervals have equal lengths L and are 
separated by a distance R > L. Then the quasiparticle model gives (see Fig. 1 for a 
plot) 


S(t) = 2Se, X 


2f, 


t < - 

L, 


— < t < — 

2 ^ ^ ^ 2 ’ 

L- 


Ft ^ ^ ^ L-\-Fi 

2 ^ h ^ 2 ’ 

L + 

1 

1 

bo|^ 

L+R ^ ^ 2L+R 

2 ^ ^ ^ 2 

L, 


t > 


(2.4) 


In particular, this result tells us that the entanglement entropy is not monotonic in time. 
We can understand the drop in the entanglement entropy in the following way. Consider 
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Figure 2. An EPR pair produced at the points marked as green at the bottom of the figure. 
When the constituent particles are at the positions marked as red at the intermediate time, 
they contribute to the entanglement entropy. At the later time when the particles are at the 
positions marked as blue, they do not contribute to the entanglement entropy. This process 
leads to a decrease in the entanglement entropy in the quasiparticle picture. 

an EPR pair that started propagating from the region between the two intervals which 
make up A at the time of quench, f = 0 (see Fig. 2). At some time t < R/2, one of 
the particles will enter region A and while the other remains in the complement. At 
that time, and for times immediately following, this pair contributes one unit to the 
entanglement entropy. However, at a later time the second particle will enter the region 
A as well, and so there is an opportunity for both particles to be inside A at the same 
time. Then they will no longer contribute to the entanglement entropy. This results in 
decreased entanglement. 

2.2 Holographic Calculation 

The time evolution of entanglement entropy after a global quench can also be studied 
using the AdS-CFT correspondence [5-7, 9, 10]. According to the AdS-CFT dictionary, 
a global quench in the boundary theory is dual to throwing a spatially homogenous and 
isotropic shell of matter into the bulk. This shell will eventually collapse to form a black 
hole, which is the gravity dual to a thermal state in the CFT. 

After the quench, the geometry of the bulk is given by the time-dependent AdS- 
Vaidya geometry (displayed as a conformal diagram in Fig. 3). As a result, the area 
of the boundary-anchored extremal surfaces, and hence the entanglement entropy of 
the boundary region, will depend on time. Though we are concerned with a (1-1-1)- 
dimensional CFT, and hence a 2-|-1-dimensional bulk, Liu and Suh were able to perform 
the calculation in arbitrary dimensions and with different types of black hole. 

The local equilibrium length is given by the horizon radius, Zh- Consider a spatial 
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Figure 3. Penrose diagram of the time dependent geometry following the quench. The red 
vertical line on the right is the AdS boundary (z = 0 in Poincare patch). The green diagonal 
line is the infalling shell, and the blue diagonal line is the horizon. The dashed curve is a late 
time extremal surface, which asymptotes to the critical surface, indicated by the solid curve. 
The linear growth of entanglement entropy comes from the portion of the extremal surface 
lying along the critical surface behind the horizon. 


interval on the boundary of length L ^ Zh- For times t ^ Zh, but smaller than L/2, 
the extremal surface in the bulk anchored to the endpoints of the interval has an area 
that grows linearly with time. Geometrically, the linear growth is tied to the existence 
of a critical extremal surface behind the black hole horizon. The extremal surface 
anchored to the boundary interval goes behind the horizon and approaches the critical 
surface. The length of the portion of the extremal surface lying along the critical surface 
increases linearly with t, which leads to a linear growth in area. At f ~ L/2, there is 
a transition (the details of which do not concern us here), after which the extremal 
surface lies outside of the horizon in the black hole portion of the geometry and is no 
longer influenced by the collapsing shell. The symmetries of this geometry ensure that 
the area is time-independent in this region, and so this represents thermal saturation 
of the entropy. 

To summarize, the vacuum subtracted entropy for a single interval as computed 
holographically is given by 


S{t) = 2Seq X 


L 
2 ’ 


t > f. 


(2.5) 


where Seq is related to the AdS-radius, Tacisj horizon radius, Zh, and Newton’s constant. 


Gn, by 


■^eq 


1 2LAdS 
AGn Zh 


( 2 . 6 ) 
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Figure 4. Here we display the extremal surfaces for two intervals. The first candidate HRT 
surface is the union of the two smaller arcs (marked in red and labeled Ai and ^ 2 )- The 


second candidate is the union of the two larger arcs (marked in green and labeled ^3 and 


^4). 


This holographic result for a single interval agrees with the predictions of the quasi¬ 
particle picture (2.3). 

We can deduce the holographic result for an arbitrary collection of intervals by 
using the answer for a single interval. For example, consider the case of a pair of in¬ 
tervals, [a;i,a; 2 ] and [ 0 : 3 ,X 4 ]. We need to find the bulk extremal surface with minimal 
area anchored to those intervals on the boundary. This is called the HRT surface [2]. 
There are two candidate HRT surfaces, which we display in Fig. 4. First, there is the 
union of the bulk extremal surfaces associated to the two intervals and [ 2 : 3 , 0 : 4 ] 

individually (^1 and A 2 in Fig. 4). But the union of the bulk extremal surfaces asso¬ 
ciated to [xi,x^ and [ 2 : 2 , 2 : 3 ] is a second choice (^3 and . 4,4 in Fig. 4). The holographic 
prescription is to compute the total area in both cases and take the minimum value. 
But in each case the extremal surfaces are just unions of extremal surfaces associated 
to intervals, and the time-dependence of those surfaces is given by (2.5). 

Let Li = X 2 — Xi, L 2 = X 4 — X 3 , and R = X 3 — X 2 . Suppose Li < L 2 . Then the hrst 
possible pair of extremal surfaces (4.i and A 2 in Fig. 4) would have a time-dependence 
given by 


^«(^)=2SeqX <f. 



(2.7) 
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while the second choice of extremal surfaces (^3 and A 4 in Fig. 4) has 


2 t, 

= 2Seq X < f + f, 

I Li-\-L2-\-2Fl 


t < — 

^ ^ ^ Li-\-L2-\-R 


( 2 . 8 ) 


At each time we take the minimum of and to get S{t). The interesting 

case is when R is the smallest of the three lengths, which means that our two disjoint 
intervals are close together. Then we have 


2 t, 


t < 


R 
2 ’ 


S{t) = 2Seq X < f + f, 


2 ' 
ni+L2 
2 ’ 


22 f —R 

o \ ^ \ o 1 


t > 


L\-\-L2—R 


In the other two cases, Li < R < L 2 and Li < L 2 < R , we have 




t < hi 


S{t) = 2Seq X<t+^, 


2 

Li -\-L/2 
2 ’ 


«>f. 


(2.9) 


( 2 . 10 ) 


A plot of the case Li = L 2 < R is in Fig. 1. Unlike the quasiparticle model, the 
holographic calculation gives a non-decreasing answer for the entropy. It is easy to see 
why this is the case. For an arbitrary boundary region, each of the candidate HRT 
surfaces is the union of a collection of extremal surfaces anchored on boundary intervals. 
But each surface anchored on a boundary interval has a non-decreasing area. Therefore 
each candidate HRT surface has a non-decreasing area, and so the true HRT surface 
has a non-decreasing area (even though the identity of the true surface may change as 
a function of time). 


3 Large-c CFT Calculation 

3.1 Vacuum State Entanglement 

In this section we review the standard machinery for calculating the entanglement 
entropy of disjoint intervals in a (l-l-l)-dimensional CFT. The idea is to use the replica 
trick [4, 14, 15] to write the entanglement entropy as a limit of correlation functions 
of twist operators. The symmetries of the CFT are used to evaluate those correlation 
functions, which lets us find the entropy. In this section, the correlation functions will 
be vacuum correlation functions, and we will review how the large-c, sparse spectrum 









assumption reproduces the holographic answer in these cases. In the following section, 
we will consider the time-dependent situation of a global quench, where the relevant 
correlation functions are those of a BCFT. 

We calculate the entanglement entropy of a region A using the replica trick, which 
realizes the entropy as a limit of traces of powers of the density matrix: 

S = lim logTrp^ (3.1) 

n-s-l 1 — n 

A standard way to compute Trp^ is via a path integral on an n-sheeted cover of the 
original surface, with branch points located at the endpoints of A. This means that the 
sheets of the cover are sewn together along A, which is the path integral representation 
of the matrix multiplication that dehnes p\. Alternatively, Aip\ can be computed as 
a certain correlation function in the theory CFT”, consisting of n copies of the original 
CFT. The twist operator Tn{x) is dehned in the CFT"^ theory as the operator which 
implements the boundary conditions of the n-sheeted cover: monodromy around the 
twist operator shifts a local operator Ok{x) in the kth copy of the CFT to the same 
operator Ok+i{x) in the /c-l-lst copy of the CFT. When the region A consists of a 
union of intervals, Trp)^ can be computed (up to a constant of proportionality) as the 
correlation function of twist operators inserted at the endpoints of the intervals. The 
twist operator Tn is inserted at all of the left endpoints, and the antitwist operator T-n 
(which sends Ok —t Ok-i) is inserted at the right endpoints. To actually evaluate these 
correlation functions, we make use of the fact that twist operators are primary with 
scaling dimension 

An = - 1). (3.2) 

where c is the central charge of the CFT. 

As an illustration, consider the region A consisting of two intervals, [a;i,a; 2 ] and 
[xz,x^, with Xi < Xj+i. Then Trp (4 can be computed as a four-point function of twist 
operators. Dehning the cross-ratio p as 

_ {X 2 - Xi){x 4 : - X 3 ) 

^ {x'i - xi){xi - X 2 y 

noting that p G [0,1], and the new coordinate tn, 

{Z- Xi){x 4 - X 3 ) 


(3.3) 


(3.4) 
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we can write the required four-point function as 


{Tn{Xi)T-n{x2)Tn{X:i)T-n{Xi)) = 


V 


(X2 - Xi)(x4 - X3) 
1 — rj 


(X 3 - X 2 )(X 4 - Xi) 


2A„ 


2Ar 


(7^(0)71„(t/)7;,(1)71„(oo)) 

(3.6) 

{Tn(0)T-„(j])T„(l)T-„(,ca)), 

(3.6) 


where the correlation function appearing on the right-hand side is defined via the limit 
{Tn{0)T-n{v)Ml)T-n{oo)) = lilU {Tn{'d)T-n{v)rn{l)T-n{w)) . (3.7) 


W^OO 

A general four-point function can be evaluated using the conformal block decomposi¬ 
tion: 


(Oi( 0)O2(^)O3(1)O4(«3)) = Cf 2 C 3 V(c./!,.{fti}.I)).F(c,/!„{/!.}. 5)). (3.8) 

P 


The sum is over all primary operators in the theory, with hp {hp) being the (anti-) 
holomorphic scaling dimension of the primary operator. This sum is sometimes called 
the s-channel decomposition, but there is also a t-channel decomposition which features 
the coefficients 0 ^ 30 ^^ instead. Evaluating this sum requires an expression for the 
conformal blocks T, which there are efficient algorithms for computing, as well as 
knowledge of the coefficients Cfj, which depend on the theory in question. We are only 
interested in correlations of twist operators, so we can set all of the hi equal to the 
same value hi = hi = h = An/2. For small r], the Taylor series of the conformal block 
gives J- = 0{t])). So at small rj, the dominant term in the sum comes from 

hp = 0, which is the identity block [16]. This is the disconnected part of the four-point 
function: keeping this term alone reduces it to a product of two two-point functions: 


{Tn{Xi)T-n{x2)Tn{x3)T-n{X4)) ^ 


1 


2 A„ 


{X2 - Xi){x4 - X3) 


= {Tn{Xi)T-n{x 2 )) ( 7 ^(Xs)?!^(0:4)) . 


(3.9) 


This is a manifestation of the cluster decomposition principle. A nontrivial fact about 
large-c CFTs with a sparse spectrum is that this is the dominant contribution even 
for finite values of rj, all the way to 77 = 1/2 [16]. This can be proved by looking at a 
large-c expansion of the conformal blocks, where it is seem explicitly that the identity 
block makes the largest contribution. Here the assumption of a sparse spectrum means 
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that the number of operators with scaling dimensions less that 0{c) should not scale 
with c. 

For T] ^ 1, the contribution of the identity block in the t-channel decomposition 

says 


{Tn{Xi)T-n{,X2)Tn{x3)T-n{x4)) ^ 


1 


2A„ 


{X3 - X2){X4- Xi) 


= {Tn{x3)T-n{x2)) {Tn{x4)T-n{,Xi)) . 


(3.10) 


Again, at large-c with a sparse spectrum this formula is expected to hold down to 
r] = 1/2. In the c —>■ cxo limit there is a sharp phase transition between the s-channel 
and t-channel results air] = 1/2. Finite-c corrections should smooth this transition, but 
consideration of those effects is beyond the scope of this work. Taking the appropriate 
limit as n —?• 1, these results show that the vacuum entanglement entropy of a pair of 
intervals at large-c is given by 


S 



{X 2 - Xi){x4 - X 3 ) 


C {X4-Xi){x3-X2) 

-log-^- 

3 


(3.11) 


where e is the UV cutoff scale. This matches the holographic answer [1]. 

When there are N intervals, we must compute a 2iV-point function of twist opera¬ 
tors. There are many possible decomposition channels, and the dramatic simplihcation 
at large-c with a sparse spectrum is that there is always some channel in which the 
identity block provides the dominant contribution, and in this channel the 2iV-point 
function decomposes as a product of N two-point functions. Taking the n —)■ 1 limit 
to extract the entanglement entropy, this precisely reproduces the Ryu-Takayanagi 
formula for arbitrary numbers of intervals [16, 17]. 


3.2 Global Quench 

Following the setup of [4, 18, 19], the global quench is effectively modeled as a BCFT 
calculation. The correlation functions we need can be computed at large-c with a sparse 
spectrum using the techniques of [16]. The leading-order holographic result is obtained 
by assuming that the only primary operator of scaling dimension less than 0(c) is the 
identity, which contributes an amount of 0(c) to the entanglement entropy. This will 
also be the dominant contribution when the spectrum is sparse, since other primary 
operators will contribute at 0(1) to the entropy. If the number of operators scales with 
c, then clearly this will compete with the identity contribution and the holographic 
result does not apply. 

We wish to compute the four-point function of the twist operators after a quantum 
quench, where we take the state to be |4/(t)) = 0(t) |\k). We are interested in times 
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well after the local equilibration time, and so it turns out to be useful to model the 
local equilibration by an initially Euclidean evolution over a small imaginary time Tq, 
which can be thought of as a regulator for the calculation. Introducing the coordinate 
z = X + ir, the path integral preparing the ket state at f = r = 0 is then a path 
integral over the strip Im^ G [—ro,0), with a boundary wavefunction at z = —iro, 
while the bra state can be obtained by integrating over the strip Imz G (0,ro]. To 
hnd the correlation functions at real times t S> tq, we would continue to path-integrate 
over a real-time contour before inserting our operators, as in the Schwinger-Keldysh 
formalism. Alternatively, we can compute correlation functions in the strip Imz G 
[—To, To] for arbitrary values of the imaginary time r and then afterward analytically 
continue the answers to real time. 

For a collection of N intervals, we need to compute a 2A^-point function of twist 
operators in the strip Imz G [—To,ro]: 


Trp”= {Tn{Xi+iT)T-n{,X 2 + iT)---Tn{x 2 N-l+ir)T-n{x 2 N + 'iT))^^^.^, (3.12) 

for arbitrary values of the imaginary time r G [—To,ro]. We will take r to be a real 
number initially, and then analytically continue r ^ it in the hnal answer to extract 
the real-time post-quench correlation functions. 

We can conformally transform the strip to the upper half plane by setting 


w{z) = exp 


TT 

2to 


(z + !To) 


(3.13) 


Then the lines z = ^iro map to the positive and negative real w axes, respectively. 
The z = ir line maps to aigw = tt /2 + t/2tq. We can write the correlation function as 

/ \ 2N 




TT 

2ro 


g2ro 


{Tn (yWi) T-n {W 2 ) ■■■Tn {w2N-l) T-n (^^^2Af))uHP ) (3T4) 


where the Wi are the images the Xj. 

The upper half plane correlation functions should be computed in the context of 
BCFT, which tells us that each primary operator in the upper half plane can be thought 
of as the product of a holomorphic operator at its location times another holomorphic 
operator at the conjugate location (reflected over the real axis) [20, 21]. Then the 2N- 
point function in the upper half-plane can be computed as a 4A^-point function in the 
full plane, which is just a vacuum correlation function. It is useful to parametrize the 
four-point function in terms of real parameters rjij, which are invariant cross-ratios 
characterizing the separation of Wi and Wj, 

WiiWjj WijWjj^ 
r]ij = 1 -^ 


(3.15) 


12 








where we have used the notation Wij = Wi — Wj and Wfj = Wi — w]. Note that rjij G [0,1]. 
Also, the UHP 2iV-point function we started with should only depend on 4iV — 3 real 
degrees of freedom (after making use of the part of the conformal symmetry which 
maps the real axis to itself), so the r]ij parameters are not all independent. They are 
still a useful parametrization, however. When r]ij ^ 1, the operators at Wi and Wj are 
much closer to their respective image points than to each other. Likewise, at r]ij ~ 0 
the operators at Wi and Wj are closer to each other than to their respective images. 
This behavior helps determine efficient OPE expansion channels, as we will see below. 

The time-dependence of the correlation function is reflected in the time-dependence 
of the r]ij cross-ratios, which in the tq —)■ 0 limit are given by 


Vijit) = 1 


2 cosh^(7rt/2ro) 

cosh(7r|a;i — Xj\/2To) -|- cosh(7rt/ro) 


1 -|- exp 



(3.16) 


We see that there are sharp transitions between rjij ^ 1 and rjij ^ 0, which occur at 
half the light-crossing time, \xi — Xj\/2. 

At early times (meaning for times less than the length scales dehning the intervals, 
but still much greater than the local equilibrium scale Tq), all of the r]ij are approx¬ 
imately equal to one. Then the upper half-plane 2iV-point function approximately 
factorizes into 2N full-plane two-point functions: this is the cluster decomposition 
limit where each operator in the upper half plane is paired with its image point in the 
lower half plane. Then we hnd 


Tip^ 


-nt/2To 

2ro 


2NA„ 


(3.17) 


This leads to a linear growth of the entanglement entropy at a rate which is N times 
as fast as it would have been for a single interval (in agreement with the quasiparticle 
picture). 

At very late times, all of the r]ij parameters are very small. One can check that 
this corresponds to a different cluster decomposition limit: now adjacent operators in 
the upper half-plane are paired with each other, and their image points are paired with 
each other. So again we hnd a product of 2N full-plane two-point functions. One can 
check that in this limit the thermal entropy formula S = SeqVh is produced. 

At intermediate times, when some of the r]ij are small and others are approximately 
equal to one, we need to be more careful. This is where the real difference between 
holographic and weakly coupled CFTs lies. To simplify the notation, we will restrict 
ourselves to iV = 1 and N = 2, though similar arguments hold for all N. 

At N = 1 there are no real surprises, but it is useful to go through it to illustrate 
some of the key points. There is only a single interval and a single cross-ratio rj, and 
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the two-point function in the upper half-plane can be written as a four-point function 
in the full plane, which we analyzed above. The conformal block decomposition of the 
four-point function implies that we can write the UHP two-point function as 

{Tn (Wi) T-n (w^2))uHP = 12A„ ' (3.18) 

Even though Fn{ri) may be a very complicated function, the cluster decomposition 
limits tell us that Th —?• 1 both when r; —)• 0 and 77 —)■ 1, since the prefactor alone 
reproduces both of those limits. This makes the global quench computation easy to 
perform when tq —)• 0 (i.e., at times and distances much greater than the equilibration 
scale), since rj makes a rapid transition from zero to one in that case. Then for arbitrary 
CFTs, the factor Fn{r]) is completely inconsequential and can be dropped; the result 
agrees with both the quasiparticle model and the holographic calculation. However, 
we would like to point out that a stronger statement can be made about large-c CFTs 
with a sparse spectrum. In the limit c —)■ cxo, we have an exact formula for Fn{ri) that 
comes from demanding that the full-plane four-point function always factorize in either 
the s-channel or f-channel: 


FnM 


T] > 1 / 2 , 
(1 — T] < 1 / 2 . 


(3.19) 


Although this formula did not make a difference for a single interval, we will now see 
how things change with multiple intervals. 

With two intervals, N = 2, the relevant correlation function is a four-point function 
in the upper half-plane. Following [22], we will parametrize the four-point function as 

(r„ {w,)T-n (u. 2 ) Tn {w,) T-n (W4))u„p = „ I ^ - |A ( 

V^12h23hl4h34/ 

(3.20) 

The function Fn (different from above, but with the same name) is in principle very 
complicated, but since we only care about rfij 1 or rjij 0, only the values of Fn at 
those points are necessary to hnd the entanglement entropy. So as time passes and the 
r]ij transition from one to zero, Fn effectively becomes a piecewise constant function of 
time. This translates into an additive term in the entropy, which is usually argued to 
be sub leading in the tq —)■ 0 limit. However, this reasoning does not apply if the leading 
constant term in Fn is zero during any given phase. Then we have to consider terms 
which are proportional to some of the nearly-vanishing rjij, which are exponentially 
small and in particular time-dependent. This is what happens for theories with a 
holographic dual. At large-c and with a sparse spectrum, the eight-point function in 
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the full plane, which is equal to the four-point function we wish to compute in the upper 
half-plane, always factorizes in some channel into a product of two-point functions [16]. 
By comparing this result to (3.20), we hnd that the effective factors are proportional 
to the rjij whenever the quasiparticle model (i.e., the prefactor in (3.20)) says that the 
entropy should decrease. In particular, the resulting time dependence of precisely 
cancels out the time dependence of the prefactor, meaning that the entropy remains 
saturated at its thermal value and does not decrease. 

4 Entanglement Tsunami 

It would be desirable to have a picture of propagating interacting quasiparticles to 
replace the free-streaming quasiparticle picture that seems to work well for free theories. 
Such a picture should reproduce the linear growth rate evident at early times, but 
avoid decreases in entropy. We have not been able to derive such a picture in terms of 
particles, but a simple heuristic which gives the correct answer is the “entanglement 
tsunami” of [9, 10], which we will elaborate on here. We will begin with a discussion 
that applies to one or two disjoint intervals, where the holographic entropy can be 
reproduced exactly, and then explain how to extend the idea to arbitrary numbers of 
disjoint intervals where we only have an upper bound on the entropy. 

4.1 One or Two Intervals 

Our picture of the entanglement tsunami is as a wave which begins at each of the 
endpoints of 4 at f = 0 and flows outward in both directions (see Fig. 5). This wave is 
not a physical wave representing the propagation of particles or energy: it is merely a 
tool for understanding the entanglement. In particular, note that the quench state is 
homogeneous, while the wave begins at particular locations picked out by the region A 
we have chosen. At any time f > 0, the wave divides the space into two regions: one 
which has already been overtaken by the wave, and one which has yet to be overtaken 
by the wave. We will call the former region the “entangled” region, and the latter region 
the “unentangled” region. This picture is reminiscent of a vacuum decay, where at f = 0 
we have tunneling events at each of the endpoints of A and the entanglement tsunami 
wavefront is like a bubble wall which converts the metastable vacuum (the unentangled 
region) into the true vacuum (the entangled region). (Though we should emphasize 
once more that there is no physical sense in which the entanglement tsunami is changing 
the vacuum; it just has the same cartoon picture.) If A has more than one endpoint, 
i.e., it is not just a half-line, then eventually the entanglement tsunami wavefronts 
starting from adjacent endpoints will collide. For all times after the collision, the entire 
interval between those two endpoints will be part of the entangled region. 
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On the left, we show the entanglement tsunami wavefront as a function of time (jagged black 
line.) The region A is marked as red. The intervals between the disconnected components of 
A are marked as blue. On the right we show the entanglement entropy as a function of time. 

We still have to give a rule for computing the entanglement entropy of the region 
A given this tsunami picture. Let the entangled region at time t be denoted by -E(t). 
When A consists of just one or two intervals, then its entanglement entropy is 

S(t) = Seq X inin(Vol(E(t) n A),Vol(E(t) n A^^)^ (4.1) 

It is not hard to see that this rule agrees with the holographic prescription for entan¬ 
glement entropy: for one interval it is trivial, and for two intervals the two options 
essentially coincide with the two possible HRT surfaces. Also, note that the entropy is 
symmetric with respect to A and A'^, which is as it should be for a globally pure state. 

Although we have emphasized that the entanglement tsunami does not represent 
the propagation of a physical excitation, there is a suggestive interpretation in which 
we do imagine a collection of excited particles living in E(t). Suppose that a hnite 
density of qubits populates the region E(t), and that those quibits are in a typical 
pure state. Then the entanglement entropy of the qubits in E(t) A A will follow the 
Page rule, meaning that their entanglement entropy will be proportional to either the 
number of qubits in E{t) fl A or E{t) n A'^, whichever is smaller [23]. This is precisely 
the entanglement tsunami prescription. 

4.2 Multiple Intervals 

When A consists of more than two intervals, there does not appear to be a simple rule 
like (4.1) which correctly reproduces the holographic answer, but right answer. One can 
attempt the following simple and natural generalization, which involves a rehnement 
of our notion of the entangled region E{t). Instead of a single region, the entangled 
region is naturally the union of disjoint intervals: 

i 
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Figure 6. Entanglement tsunami for many intervals. The region A is marked as red. The 
intervals between the disconnected components of A are marked as blue. Note that at time 
ti, E{t) consists of four disconnected components (orange solid lines) separated by the entan¬ 
glement tsunami wavefront (jagged black line). But at time t 2 , the first pair and second pair 
have merged, leaving two disconnected regions. At time there is only a single connected 
region. 

where each Ei{t) is an interval representing a single connected component of E{t). 
The number of Ei{t) regions changes with time. For N intervals, E{t) has iV + 1 
connected components at early times (each being a small neighborhood around an 
endpoint of an interval) and only a single connected component at late times. A collision 
of entanglement tsunami wavefronts indicates that two of the Ei{t) are merging, and 
thereafter will be treated as a single unit. We illustrate this behavior in Fig. 6. 

With this rehnement of the entangled region, we have the following upper bound 
on the holographic entanglement entropy^: 

S{t) < Seq^min(^Vol(Ei(t) n A),Vol(Ei(t) n A")). (4.3) 

i 

In other words, we have a separate minimization problem for each connected component 
of the entangled region, and at the end we add them all up. This rule is numerically 
the same as (4.1) for one or two intervals, and so there is equality. For more than two 
intervals, (4.3) corresponds to the area of one of the candidate HRT surfaces, though 
not necessarily the minimal one. Therefore it represents only an upper bound on the 
entropy. Forthcoming work from Casini, Liu, and Mezei discusses such an upper bound 
beyond the context of holography [24]. We also note that (4.3) is manifestly symmetric 
with respect to A and A'^, even though each minimization subproblem is free to choose 
to use region A or independently for its contribution to the entropy. 

There is also still a suggestive qubit picture for the right-hand side of (4.3). Now 
instead of one collection of qubits in E{t), there is an independent collection of qubits 
in each of the F^i(f), and the total qubit state is a product state over the connected 

^We would like to thank Mark Mezei for providing us with an example of three disjoint intervals 
where equality fails. 
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components. We can imagine that the qubits are somehow being emitted at the end¬ 
points of the intervals, and when two qubit chains come into contact (i.e., when two 
wavefronts merge), the two qubit chains undergo a very rapid mixing and end up in a 
typical pure state of the composite system. Even though this picture only provides an 
upper bound on the entropy, it could be that a similar picture is accurate for the real 
system. 

5 Discussion 

5.1 Summary 

We have seen that for a (l-l-l)-dimensional CFT in the limit of large-c, with a sparse 
spectrum of operators, the time-dependence for the entanglement entropy of multiple 
intervals is not correctly captured by the quasiparticle model of [3, 4]. The correct 
answer, which is correctly reproduced holographically, is non-decreasing with time, 
whereas the quasiparticle model features crops in the entropy. A heuristic model which 
gives the correct time evolution is provided by the entanglement tsunami, a picture 
originally introduced in [9, 10] which we have developed into a rule for calculating the 
entanglement. There are many unanswered questions about this model, which we will 
now discuss. 

5.2 Higher Dimensions 

In higher dimensions, as in (l-l-l)-dimensions, the entanglement entropy for disjoint 
regions in the quasiparticle picture will experience periods of decrease that are not 
present in the holographic calculation. The entanglement tsunami picture will solve 
this problem since, like the holographic calculation, it does not allow for decreases in 
the entropy. There remains the puzzle of the velocity of the entanglement tsunami 
wavefront. In (l-l-l)-dimensions, it seems natural to say that the wavefront moves 
at the speed of light. But this would cause the entanglement to grow more quickly 
than the holographic calculation indicates. At the same time, a wavefront made up 
of quasiparticles traveling in random directions moves too slowly. A held-theoretic 
derivation of the entanglement tsunami velocity would do much to clarify the picture. 

5.3 Interacting Quasiparticles? 

A natural guess for hxing up the quasiparticle picture is simply to add interactions. 
Instead of free-streaming, the particles would bump into each other and generate multi¬ 
partite entanglement. It is not clear whether a simple picture of this type can accurately 
reproduce the holographic calculation, but the entanglement tsunami prescription may 
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provide a useful starting point. In the entanglement tsunami, it appears as though 
there are clouds of particle being emitted by the endpoints of the the region A (or, in 
higher dimensions, by the boundary of A), beginning at the quench time. The unphys¬ 
ical aspect of this model is that it is not homogeneous: the boundaries of A are not 
preferred locations in the state, so it doesn’t make sense to say that they are the sources 
of particles. Instead, the state should describe a hnite density of interacting particles 
everywhere at once. However, since we are asking about the entanglement entropy of A, 
there may be a sense in which the particles in a neighborhood of the boundary of A are 
the only ones that matter. Then the entanglement tsunami would represent a growing 
sphere of influence, picking out the parts of the homogenous matter distribution which 
determine the entanglement. We leave an exploration of this possibility to future work. 

Acknowledgments 

We would like to thank C. Akers, R. Bousso, T. Faulkner, Z. Fisher, T. Grover, 
V. Hubeny, A. Kitaev, J. Koeller, M. Mezei, R. Myers, S. Shenker, M. Smolkin, D. Stan¬ 
ford, J. Suh, and C. Zukowski for useful discussions and correspondence, as well as all of 
the participants in the “Quantum Gravity Foundations: UV to IR” and “Entanglement 
in Strongly-Gorrelated Quantum Matter” programs at KITP, Santa Barbara. This work 
was supported in part by the Berkeley Genter for Theoretical Physics, by the National 
Science Foundation (award numbers 1214644 and 1316783), by fqxi grant RFP3-1323, 
and by the US Department of Energy under Gontract DE-AG02-05GH11231. 

References 

[1] S. Ryu and T. Takayanagi, “Holographic derivation of entanglement entropy from 
AdS/CFT,” Phys.Rev.Lett. 96 (2006) 181602, arXiv:hep-th/0603001 [hep-th] . 

[2] V. E. Hubeny, M. Rangamani, and T. Takayanagi, “A Covariant holographic 
entanglement entropy proposal,” JHEP 0707 (2007) 062, arXiv:0705.0016 [hep-th]. 

[3] P. Calabrese and J. L. Cardy, “Evolution of entanglement entropy in one-dimensional 
systems,” J.Stat.Mech. 0504 (2005) P04010, arXiv: cond-mat/0503393 [cond-mat] . 

[4] P. Calabrese and J. Cardy, “Entanglement entropy and conformal field theory,” 

J.Phys. A42 (2009) 504005, arXiv: 0905.4013 [cond-mat. stat-mech] . 

[5] J. Abajo-Arrastia, J. Aparicio, and E. Lopez, “Holographic Evolution of Entanglement 
Entropy,” JHEP 1011 (2010) 149, arXiv: 1006.4090 [hep-th]. 


19 


[6] V. Balasubramanian, A. Bernamonti, J. de Boer, N. Copland, B. Craps, et al, 
“Thermalization of Strongly Coupled Field Theories,” 

Phys.Rev.Lett. 106 (2011) 191601, arXiv: 1012.4753 [hep-th], 

[7] V. Balasubramanian, A. Bernamonti, J. de Boer, N. Copland, B. Craps, et al., 
“Holographic Thermalization,” Phys.Rev. D84 (2011) 026010, 

arXiv:1103.2683 [hep-th]. 

[8] T. Hartman and J. Maldacena, “Time Evolution of Entanglement Entropy from Black 
Hole Interiors,” JHEP 1305 (2013) 014, arXiv : 1303.1080 [hep-th]. 

[9] H. Liu and S. J. Suh, “Entanglement Tsunami: Universal Scaling in Holographic 
Thermalization,” Phys.Rev.Lett. 112 (2014) 011601, arXiv : 1305.7244 [hep-th]. 

[10] H. Liu and S. J. Suh, “Entanglement growth during thermalization in holographic 
systems,” Phys.Rev. D89 (2014) no. 6, 066012, arXiv: 1311.1200 [hep-th]. 

[11] V. Balasubramanian, A. Bernamonti, N. Copland, B. Craps, and F. Galli, 
“Thermalization of mutual and tripartite information in strongly coupled two 
dimensional conformal field theories,” Phys.Rev. D84 (2011) 105017, 
arXiv:1110.0488 [hep-th]. 

[12] C. T. Asplund and A. Bernamonti, “Mutual information after a local quench in 
conformal field theory,” Phys.Rev. D89 (2014) no. 6, 066015, 

arXiv:1311.4173 [hep-th]. 

[13] A. Allais and E. Tonni, “Holographic evolution of the mutual information,” 

JHEP 1201 (2012) 102, arXiv: 1110.1607 [hep-th]. 

[14] P. Calabrese and J. L. Cardy, “Entanglement entropy and quantum field theory,” 
J.Stat.Mech. 0406 (2004) P06002, arXiv:hep-th/0405152 [hep-th]. 

[15] P. Calabrese, J. Cardy, and E. Tonni, “Entanglement entropy of two disjoint intervals 
in conformal field theory,” J.Stat.Mech. 0911 (2009) PllOOl, 

arXiv:0905.2069 [hep-th]. 

[16] T. Hartman, “Entanglement Entropy at Large Central Charge,” 
arXiv:1303.6955 [hep-th]. 

[17] T. Faulkner, “The Entanglement Renyi Entropies of Disjoint Intervals in AdS/CET,” 
arXiv:1303.7221 [hep-th]. 

[18] P. Calabrese and J. Cardy, “Quantum Quenches in Extended Systems,” 

J.Stat.Mech. 0706 (2007) P06008, arXiv: 0704.1880 [cond-mat. stat-mech] . 

[19] P. Calabrese and J. L. Cardy, “Time-dependence of correlation functions following a 
quantum quench,” Phys.Rev.Lett. 96 (2006) 136801, 

arXiv:cond-mat/0601225 [cond-mat] . 


20 


[20] J. L. Cardy, “Boundary conformal field theory,” arXiv:hep-th/0411189 [hep-th] . 

[21] P. Di Francesco, P. Mathieu, and D. Senechal, “Conformal field theory,”. 

[22] A. Coser, E. Tonni, and P. Calabrese, “Entanglement negativity after a global 
quantum quench,” J.Stat.Mech. 1412 (2014) no. 12, P12017, 

arXiv:1410.0900 [cond-mat.stat-mech] . 

[23] D. N. Page, “Average entropy of a subsystem,” Phys.Rev.Lett. 71 (1993) 1291-1294, 
arXiv:gr-qc/9305007 [gr-qc] . 

[24] H. Casini, H. Liu, and M. Mezei, “To appear,”. 


21 


